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ABSTRACT 

This paper investigates the effect of differential aperture loss with SDSS fibers 
and examines whether such selection bias would result in the observed correlation 
between rest-frame absorption equivalent width of Mgll absorbers, Wr{2796), and 
mean associated [O II] luminosity, L[o n] , in SDSS QSO spectra. We demonstrate based 
on a Monte Carlo simulation that the observed Wr{2796) vs. i^oii] correlation of 
Mg II absorbers can be well-reproduced, if all galaxies found in deep surveys possess 
extended Mg II halos and if the extent of Mg II halos scales proportionally with galaxy 
mass as shown in previous studies. The observed correlation can be explained by a 
combination of (1) the known Wr(2796) vs. p anti-correlation in galaxy and Mgll 
absorber pairs and (2) an increasing aperture loss in the 3" diameter SDSS fiber for 
galaxies at larger p. Galaxies at larger projected distances produce on average weaker 
Mgll absorbers and weaker (or zero) ipii] in SDSS QSO spectra. We show that such 
correlation diminishes when larger fibers are adopted and is therefore not physical. 
While under a simple halo model the majority of Mgll absorbers do not directly 
probe star- forming disks, they trace photo- ionized halo gas associated with galaxies. 
We show that because of the scaling relation between extended gas cross-section and 
galaxy mass, the number density evolution of the Mgll absorber population as a whole 
provides a good measure of the cosmic star formation history. 
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1 INTRODUCTION 

Absorption line spectroscopy is a powerful tool for study- 
ing the structure of the distant universe. By observing the 
absorption features imprinted in the spectra of background 
QSOs, we can study otherwise invisible gaseous structures 
to redshift as high as background QSOs can be. The Sloan 
Digital Sky Survey (SDSS; York et al. 2000) has yielded 
optical spectra of ^ 106 k QSOs at z — 0.065 — 5.46 (Schnei- 
der et al. 2010). This unprecedentedly large sample of QSO 
spectra offers a rich resource for studying the distant uni- 
verse using absorption spectroscopy. For example, MgllAA 
2796,2803 absorption doublets are commonly seen in QSO 
spectra obtained using ground-based spectrographs. They 
provide a uniform probe of intervening gas over a broad 
redshift range from z = 0.4 to z « 2.5. Several groups have 
conducted systematic searches of Mgll absorption features 
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in SDSS QSO spectra, producing a large sample of these 
absorbers for constraining their statistical properties (e.g. 
Bouche et al. 2004; Nestor et al. 2005; Prochter et al. 2006; 
York et al. 2006; Quider et al. 2011). 

Mgll absorbers are routinely found in the vicinities 
( ^ 100 kpc) of distant galaxies and provide a convenient 
probe of galactic halos at high redshifts (e.g. Bergeron 1986; 
Steidel et al. 1994; Kacprzak et al. 2008; Chen & Tinker 
2008; Gauthier & Chen 2011). However, their physical ori- 
gin, whether the absorbers arise in infalling clouds, outflows 
from nearby star-forming regions, or a combination thereof, 
is unclear. The ubiquitous presence of blueshifted Mg II ab- 
sorption features along the sightlines into the star-forming 
regions of z ~ 1 galaxies indicates that starburst driven out- 
flows are common in distant galaxies and that outflows con- 
tribute to some fraction of Mg II absorbers uncovered along 
random QSO sightlines (e.g. Weiner et al. 2009; Rubin et al. 
2010). Although the extent of galactic-scale winds around 
these galaxies is unknown, the observed self-absorption of 



2 Lopez & Chen 



Mgll doublets in distant star-forming galaxies has moti- 
vated recent works that attribute all strong Mg II absorbers 
(rest-frame absorption equivalent width Wr (2796) > 0.5 A) 
along QSO sightlines to outflows (e.g. Menard & Chelouche 
2009; Chelouche & Bowen 2010; Menard et al. 2011). 

Using a sample of 8500 Mgll absorbers at z = 0.4 — 1.4 
from SDSS DR4 quasar spectra, Menard et al. (2011; here- 
after Mil) observed a strong correlation between VKr(2796) 
and their associated median [O II] luminosity per unit area, 
^^^^ ^, in stacked QSO spectra. The observed Wr{2796) 

vs. y^, correlation is characterzed by 



^(O II] 



W,(2796) 



lA 



(1) 



where a = (1.48 ± 0.18) x lO^'^ ergs"^ kpc^ and 6 = 1.75 ± 
0.11. This observed correlation applies to Mgll absorbers of 
iyr(2796) — 0.7 — 6 A. A similar trend has also been men- 
tioned in Noterdaeme et al. (2010), but these authors did 
not find a correlation between 1^^(2796) and I/[o ii] in [O II]- 
emission selected Mgll absorbers. In addition, Noterdaeme 
et al. (2010) pointed out that part of their strong Mg II ab- 
sorbers arise in low -£/[oii] galaxies. Because [O II] luminosity 
L[oii] provides a measure of current star formation rate (e.g. 
Kennicutt 1998), Mil interpreted the observed strong corre- 
lation as Mg II absorbers tracing on-going star formation. In 
addition. Mil showed that the frequency distribution func- 
tion of Mgll absorbers and the [Oil] luminosity function 
share similar shape and amplitude, and that the number 
density evolution of Mgll absorbers follows the cosmic star 
formation history. Combining these empirical correlations, 
the authors argue that outflows are the mechanism respon- 
sible for the observed Mg II absorption in QSO spectra. 

The conclusion of outflows being responsible for the ob- 
served Mgll absorbers appears to be discrepant from pre- 
vious findings that VKr(2796) correlates more strongly with 
galaxy mass and weakly with galaxy colors or recent star for- 
mation history (e.g. Steidel et al. 1994; Chen et al. 2010a,b). 
Such conclusion also makes it difficult to understand the 
identifications of strong Mgll absorbers in the vicinities of 
quiescent galaxies (e.g. Gauthier et al. 2010; Bowen & Che- 
louche 2011; Gauthier & Chen 2011). 

While Mil presented a clever approach to estimate the 
co-moving star formation rate density p* based on the ob- 
servations of Mg II absorbers in SDSS QSO spectra, it is im- 
portant to understand the underlying factors that shape the 
observed strong Wr{2796) vs. X/Ljq j correlation in SDSS 
data. We note that the SDSS fibers have a finite size of 3" 
diameter on the sky, which corresponds to projected phys- 
ical separations of 11 — 18 kpc at z = 0.4 — 1.4. The 
observed W^r(2796) vs. p anti-correlation (e.g. Chen et al. 
2010a) implies that galaxies at larger projected distances 
produce on average weaker Mgll absorbers and lower (or 
zero) i[oii] (when the star-forming disks occur at angular 
distances 9 > 1.5", or p > 6 - 9 /i"^ kpc, of the QSO 
sightlines; see Figure 1). This selection bias strenthens the 
apparent correlation of Equation (lj3- A. similar point on 



^ Menard et al. (2011) applied a fiber selection correc- 
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Figure 1. Cartoon illustrating the fiber selection bias of interven- 
ing galaxies in SDSS QSO spectra. Intervening galaxies at difl'er- 
ent projected distances to the QSO sightline may imprint absorp- 
tion and/or emission features in the spectrum of the background 
QSO. Galaxies with optical disks overlapping a SDSS QSO fiber 
are expected to imprint [O II] emission features in the QSO spec- 
trum. Galaxies with gaseous halos intercepting the QSO sightline 
are expected to imprint Mgll absorption features in the QSO 
spectrum. In this picture, the optical disks of galaxies C and D 
are covered by the fiber (in part or in full), galaxies A and B are 
not. Only C and D are expected to contribute to the [Oil] emis- 
sion line seen in the stacked QSO spectrum. On the other hand, 
the halos of galaxies B, C, and D all intercept the QSO sightline, 
while A does not. Only B, C, and D are therefore expected to con- 
tribute to the Mgll absorption statistics. But because the QSO 
sightline impacts the outer part of galactic halo B, the pathlength 
through the halo is shorter and the Mgll absorbing strength is 
expected to be weaker (e.g. Chen et al. 2010a). 



the potential missed galaxy light in the SDSS fibers has also 
been made in Noterdaeme et al. (2010). 

To investigate the selection bias as a result of differential 
aperture loss, we have carried out a Monte Carlo simulation 
study. Adopting a empirical model for describing gaseous 
clouds around galaxies from Chen et al. (2010a) and empir- 
ical relations for describing the luminosity and size distribu- 
tions of the general galaxy population, we demonstrate in 
this paper that the observed Wr(2796) vs. -L[oii] correlation 
can be well re-produced without any fine-tuning after ac- 
counting for the differential aperture loss of galaxy fluxes in 
the SDSS fiber. While under the simple halo model the ma- 
jority of Mgll absorbers do not directly probe star-forming 
disks like high column density damped Lya absorbers (e.g. 
Wild et al. 2007), they do probe photo-ionized clouds around 
distant galaxie^j If more massive galaxies are surrounded by 
more extended Mgll absorbing gas, then the number den- 
sity evolution of Mgll absorbers naturally follows the cosmic 
star formation history. Because models which do not require 



^[OII] (2)/^flbcr(2), where Afn^g,.(z) represents the corresponding 
physical area of SDSS fiber at redshift z. However, as illustrated 
in Figure 1, the reduction from total observed fluxes to fiux sur- 
face density does not correct for the differential aperture loss of 
galaxy fluxes in SDSS QSO spectra. 

^ Note that the empirical 14^(2796) vs. p relation of Chen et al. 
(2010a) was established based on a sample of galaxies at close 
projected distances to a QSO sightline, including those that may 
produce a damped Lya absorption feature in the QSO spectrum. 
We therefore expect that absorber samples generated based on the 
mean Wr{2796)-p relation and the observed scatter include con- 
tributions from star-forming disks. Based on the observed number 
densities of damped Lyo absorbers (Rao et al. 2006) and Mgll 
absorbers of VKr(2796) > 0.5 A (Nestor et al. 2005) at ^ < 1, we 
estimate that no more than 20% of these strong Mgll absorbers 
have contributions from star-forming disks 
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outflows can also reproduce the empirical correlations be- 
tween absorber abundances and star formation properties, 
we caution drawing conclusions in favor of an outflow ori- 
gin for QSO absorbers based on these simple correlations. 
Throughout this paper, we adopt a ACDM cosmology with 
Qm = 0.3 and Qa = 0.7 and a dimensionless Hubble con- 
stant of ft = i?o/(100kms"^ Mpc"^). 



2 MOCK GALAXY AND ABSORBER 
CATALOGS 

To sirrmlate the SDSS observations, wc flrst generate a mock 
catalog of 200,000 gala^xies distributed uniformly within 
« 10" angular radius of a QSO sightline at redshift between 
z = 0.4 and z = 1.4. The large number is necessary to pro- 
vide a representative sampling of a wide range of galaxy 
properties (such as size and luminosity) and to minimize 
statistical noise in our simulations. The redshift range is se- 
lected to match the study of Mil. Over the redshift range of 
z — 0.4— 1.4, the angular radius of 10" corresponds roughly 
to p = 37 - 59 kpc. 

The working hypothesis of this exercise is that all galax- 
ies are surrounded by extended gaseous halos. This hypothe- 
sis is supported by empirical observations that reveal a high 
gas covering fraction around galaxies of a wide range of lu- 
minosity and color (e.g. Chen et al. 2010a). The gaseous 
halos are expected to produce Mgll absorption features in 
the spectrum of a background QSO when intercepting the 
QSO sightline. For every Mgll absorber in the mock cat- 
alog, the location and intrinsic properties of the absorbing 
galaxy are known. It is therefore possible to compute the 
emission fluxes of the absorbing galaxies recorded in QSO 
spectra. 

2.1 The Mock Galaxy Catalog 

We first generate random galaxies following a Schechter 
probability function of their rest-frame absolute B-band 

magnitude Mb, 

P(Mb) oc 10 ^ A -r ; exp(— 10 ^ '), (2) 

where a is the faint-end slope of the galaxy luminosity func- 
tion. We adopt a = —1.3 for the faint-end slope (e.g. Faber 

et al. 2007) and Mb, - 5 log ft = -19.8 following Chen 
et al. (2010a). The mock galaxy sample spans a lumionsity 
range between 0.01 L, and 10 L,. As discusse below, we also 
repeat the Monte Carlo simulations with different faint-end 
slop values, a = —1.2 and a = —1.4 and our findings remain 
the same. 

Next, we determine the optical size of each galaxy in 
the mock catalog according to the luminosity-size relation 
from Cameron & Driver (2007), 

log (7?5o) = -0.1 Mb - 1.35 (3) 

where Rr,Q is the half-light radius, the radius within which 
the galaxy emits half of its total flux. For a galaxy of A'Ib, 
the half-light radius is drawn from a Gaussian distribution 
about the mean luminosity-size relation with a 1-a width of 
0.4 dex in the log R50 space. To include the majority of the 
galaxy light, the radius that contains 90% of the total flux 



i?90 is more relevant in our study. We therefore convert R50 
to -R90, assuming an exponential surface brightness profile. 

Next, we determine the [OH] luminosity I/[oii] for each 
galaxy assuming that the rest-frame -B-band fluxes trace 
the mean profile of [Oil] emission. The expected L[oii] is 
then calculated according to the correlation between Mb 
and L[oii] found in deep survey data. Based on the study 
of Zhu et al. (2007), we find that the Mb vs. log L[oii] 
correlation can be characterized by 

log [L[oii]/(ergs-^)] = 35 - 0.3 Ms (4) 

with a 1-a scatter of aiog i[oii] ~ ^-^ Note that Equa- 
tion (4) is an empirical relation between observed quantities. 
It does not include dust extinction corrections. The inferred 
L[oii] can therefore be directly compared to what is ob- 
served in the stacked SDSS QSO spectra. For a galaxy of 
Mb, L[ou] is drawn from a Gaussian distribution about the 
mean A/s-I/[oii] relation with a 1-a width of 0.3 dex in the 
log I/[oii] space. 

To examine whether or not our mock galaxy sample is 
representative of the field galaxy population, we present the 
[O II] luminosity distribution of our mock galaxy sample in 
Figure 2 along with the observed [O II] luminosity function 
of 2; ^ 1 galaxies from different deep surveys (Takahashi et 
al. 2007; Zhu et al. 2009). All the empirical measurements 
have been converted to have the same cosmological parame- 
ters adopted in our analysis. Figure 2 demonstrates that we 
reproduce the observed [O II] luminosity function with the 
mock galaxy sample for the luminous galaxy population. 
At the faint-end, the observations suffer from survey incom- 
plctness (e.g. Zhu et al. 2009) and therefore the observed 
space density represents a lower limit to the underlying faint 
galaxy population. 

Next, we model each galaxy in the mock catalog as a 
round disk and randomly place the galaxy within « 10" of a 
QSO sightline and at redshift between z — 0.4 and z — 1.4. 
We compute the corresponding physical projected distance 
of the galaxy based on the redshift and angular distance to 
the QSO. For each disk, we also assign a random inclination 
angle with respect to the observer and a random position 
angle of the major axis of the inclined disk with respect to 
the line connecting the galaxy and the QSO. 

To determine the fraction of the luminous disk that 
overlaps the SDSS fiber centered at the QSO, we adopt the 
size of each disk galaxy Ego and its relative distance and 
orientation to the QSO. Wc de-project the inclined disk and 
determine the galactocentric distance R of each point (x, y) 
in the disk according to 

R'^{x,y) = r^{x,y) [l -|- cos^ (0 -|- a) tan^ i] (5) 

where r is the projected radius on the plane of the sky, a is 
the position angle of the major axis of the inclined disk, 9 
is the azimuthal angle of point (a;, y) from the major axis, 
and i is the inclination angle of the disk. The expected [O II] 
emission from the galaxy in the QSO spectral data is then 
computed by integrating all the light within Rqo of the disk 
that falls in the 3" diameter fiber. We consider two different 
surface brightness profiles: (1) a fiat, top- hat profile and (2) 
an exponential profile for this calculation. If the luminous 
disk does not overlap the 3" diameter fiber, then we set 
^[oii] = 0. Finally, we divide the computed i[oii] by the 
physical area of a 3" diameter fiber at the redshift of the 
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Figure 2. The [Oil] luminosity distribution of the mock galaxy 
sample described in Section 2. The histograms display the sim- 
ulated data for different faint-end slope of a = —1.2 (dotted 
histogram) a = —1.3 (solid histogram), and a = —1.4 (dashed 
histogram) . The [Oil] luminosity function measured for 2 ~ 1 
galaxies in DEEP2 (closed circles; Zhu et al. 2009) and for z ~ 1.2 
galaxies in the COSMOS field (closed squares; Takahashi et al. 
2007) and the SDF (closed triangles; Takahashi et al. 2007) are 
also shown for comparison. 

mock galaxy in order to calculate its [O II] luminosity surface 
density , that would be recorded in an SDSS QSO 
fiber. 

2.2 The Mock Mg II Absorber Catalog 

The mock Mgll catalog is formed by calculating the ex- 
pected Mgll absorption strength in the spectrum of the 
background QSO for every galaxy in the mock galaxy cata- 
log. To determine the associated Mgll absorption strength 
of a galaxy in the QSO spectrum, we adopt the uniform 
gaseous halo model of Chen et al. (2010a). Under this model, 
the Mgll absorber strength for a galaxy at projected dis- 
tance p is characterized by a mean relation of 

1^.(2796) = tan-i J f (6) 

and dispersion criog = 0.233, where -Rgas is the gaseous 
radius of the Mgll halo, an is the core radius and is an ~ 
0.2 Rgas and log Wo = 1.24 ± 0.03. The gaseous radius is 
determined based on the galaxy B-band luminosity following 
a power-law model, 



Rgiis, \ Lb, J 

for which Chen et al. (2010a) found a best-fit character- 
istic radius of -Rgas, = 75 kpc for an L, galaxy and 



a scaling index of j3 = 0.35. The scaling relation between 
gaseous radius and galaxy B-band luminosity is understood 
as more massive galaxies are surrounded by more extended 
halos (e.g. Tinker & Chen 2008; Chen et al. 2010b). In ad- 
dition to the spatial profile of the extended Mg II gas, Chen 
et al. (2010a) also measured a high gas covering fraction k 
within Rgas. Specifically, they found a mean covering frac- 
tion of Ko.3 ~ 100% for absorbers of W^r(2796) ^ 0.3 A at 
p ^ 0.4 7?gas and kq.s ~ 70% at p < -Rgas- The gas covering 
fraction increases with decreasing W^r(2796) threshold and 
decreases with p (see Figure 10 in Chen et al. 2010a). 

For each simulated galaxy, we compute the expected 
Wr (2796) by randomly sampling the mean relation Equa- 
tions (6) within the observed scatter criog w,, . For simulated 
dwarf galaxies with -Rgas < p, we set Wr (2796) — 0. In- 
cluding the large scatter (criog vi/,, ) in the computation of 
VKr(2796) allows the possibility of galaxies at small p pro- 
ducing Mg II absorbers that are much weaker than the mean. 
Although every galaxy at p < -Rgas in the mock sample is 
expected to produce a Mg II absorber, the gas covering frac- 
tion measured at a given Wr{2796) threshold is not 100%. 

For a mock sample of 200,000 galaxies, the procedure 
described above produces a mock sample of ~ 23, 000 Mg II 
absorbers of H/r(2796) ^ 0.3 A at z = 0.4 - 1.4. The fre- 
quency distribution function of the mock Mgll absorber 
sample is presented in Figure 3. For comparison, we also 
include in Figure 3 the model from Prochter et al. (2006) 
that best describe the incidence of Mgll absorbers iden- 
tified a,t z — 0.4 — 2.3 in their blind survey. Aside from 
the offset in the normalization. Figure 3 shows that the 
frequency distribution function of the mock Mg II absorber 
sample follows the shape of the best-fit model of Prochter 
et al. (2006). The difference in the incidence of mock Mgll 
absorbers (up to 20%) for different adopted faint-end slope 
of the galaxy luminosity function is understood by the rela- 
tively steep dependence of halo gas cross section on galaxy 
luminosity, cTgas = ^r^gas oc L'^^ from Equation (7). Adopt- 
ing known galaxy luminosity functions at 2: < 1 from Faber 
et al. (2007) and the scaling relation of Equation (7), we 
expect to find (following Equation 8 below) a number den- 
sity of n{z) ~ 0.8 (0.2) per line of sight for Mgll absorbers 
of Wr(2796) > 0.3 (1.0) A. Our model therefore well repro- 
duces the observe absorber statistics at 2: < 1 (e.g. Nestor 
et al. 2005; Prochter et al. 2006). 



3 ANALYSIS 

Using the mock galaxy and Mg II absorber catalogs, we pro- 
ceed to examine the relation between Wr(2796) and -L[oii] 
of Mgll absorbers. For each Mgll absorber in the mock cat- 
alog, we know from our simulation the location and luminos- 
ity of the absorbing galaxy. We have also calculated accord- 
ing to the procedures described in Section 2.1 the fraction of 
galaxy light that would be recorded in the SDSS QSO fibers. 
We can therefore directly compare our simulation data with 
the observations of Mil. 

We present in Figure 4 the distribution of ^,-(2796) 
and X/-t[o J mock Mgll absorber sample. The top 

panel shows the Wr(2796)-y^ j distribution of individ- 
ual Mgll absorbers, including detections (dots) and non- 
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Figure 3. The frequency distribution function of tlie mock Mgll 
absorber sample. As in Figure 1, the dotted, solid, and dashed his- 
tograms represent mock absorber samples from simulated galaxies 
of faint-end slopes of a = —1.2, —1.3, and —1.4, respectively. For 
comparison, we also include the observed absorber frequency dis- 
tribution function from Prochter et al. (2006) for Mgll absorbers 
found at 2 = 0.4 - 2.3. 



detections (indicated by the arrows). The bottom panel 
shows the fraction of Mg II absorbers for which the absorbing 
galaxies are not expected to overlap the fiber and therefore 
have , = in the QSO spectra. It is clear that a srow- 

ing fraction of galaxy light is missed in the QSO spectra for 
absorbers of decreasing strength due to the Wr(2796) vs. p 
anti-correlation (e.g. Chen et al. 2010). 

To reproduce the observations in SDSS QSO spectra, 
we divide the mock catalog of 23,000 Mgll absorbers into 
subsamples according to their absorption strengths. Follow- 
ing the proceduce described in Mil, who formed a median 
QSO spectrum at the rest-frame of the Mgll absorbers and 
measured the associated [Oil] emission line flux, we com- 
pute the median value of [Oil] luminosity surface density 
j)med for all Mgll absorbers in each Wr{2796) bin, 

including those with X/-^[oii] = 0- To estimate the scatter, 
we repeat the Monte Carlo simulation 100 times to generate 
100 mock samples of galaxies and Mgll absorbers. We mea- 
sure the l-CT dispersion in (X^^^ ^ )mcd among the 100 mock 
samples. The stars and the associated errorbars in Figure 
4 represent the computed (X^^ )mcd and the associated 

dispersion in each VKr(2796) bin. 

For comparison, the best-fit power-law model of the 
observed W, (2796)-(^^^ j)mcd correlation (Equation 1) 
from Mil is included in Figure 4 as the dash-dotted line. 
The Mil model is found to match well with the simulated 
data. However, both fall below the locus defined by individ- 
ual dots because of non-detections quantified in the bottom 



Figure 4. The Wr(2796) vs. Y^, distribution of the mock 

galaxy and absorber samples, in comparison to observations of 
Mil. The dots in the top panel represent individual Mgll ab- 
sorbers with detectable emission fluxes in the QSO spectra from 
their absorbing galaxies. Recall that these are calculated by di- 
viding the expected [Oil] line flux in the fiber with the corre- 
sponding physical area of the 3" fiber at the redshift of the ab- 
sorber. The arrows at the bottom row indicate the presence of 
absorbing galaxies with ~ ^ (non-detections), because 

their optical disks do not overlap with the fiber. For complete- 
ness, the bottom panel shows the fraction of Mg II absorbers that 
have ^2li j ~ ^' '^^^ (red) stars in the top panel represent 
the median values of ^2l^ j '^^^ Mgll absorbers (including 
detections and non-detections) in each iyr(2796) bin with the 
l-cr dispersion indicated by the errorbars. The dash-dotted line 
shows the best-fit power-law model (Equation 1) of the observed 
Wr(2796)-(^^^^^^^>nicd correlation from Mil. The solid (blue) 

circles represent the mean ^ averaged over all Mgll ab- 

sorbers in each Wr(2796) bin. Note that the stars fall below the 
locus defined by individual dots because of non-detections quan- 
tified in the top panel. 

panel. We also calculate the mean [Oil] flux X/l^q j 
pected in the stacked SDSS QSO spectra, averaged over all 
Mgll absorbers in each bin including non-detections. The 
results are shown in solid circles in Figure 4. 



4 DISCUSSION AND SUMMARY 

We have carried out a Monte Carlo study to investigate the 
effect of differential aperture loss of extended emission from 
intervening galaxies in SDSS flbers. We generate a mock 
galaxy sample based on known empirical relations that de- 
scribe the luminosity and size distribution of the general 
galaxy population uncovered in deep surveys. The mock 
galaxy sample is accompanied by a mock Mgll absorber 
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sample that is generated based on a simple assumption that 
extended gaseous halos are a common and generic feature 
of distant galaxies and the gaseous extent scales propor- 
tionally with galaxy mass. For each absorber in the mock 
sample, the luminosity and projected distance of the absorb- 
ing galaxy are known, allowing us to make predictions for 
the observed relation between Wr-(2796) and the associated 
L[oii], in SDSS QSO spectra. 

Our study shows that combining the known VFr(2796) 
vs. p anti-correlation of Mg Il-galaxy pairs (e.g. Chen et 
al. 2010a) and differential fiber selection of Mgll absorbing 
galaxies (Figure 1) reproduces the observed Wr{2796) vs. 
L[oii] correlation in the SDSS QSO spectra without addi- 
tional fine-tuning or scaling. The results of our study in- 
dicate that the observed 14^^(2796) vs. Z/[oii] correlation of 
Mgll absorbers in SDSS data is likely due to a differential 
fiber selection bias and does not provide a physical under- 
standing of the origin of the Mgll absorber population. 

On the basis of the the observed Wr(2796) vs. -£/[oii] 
correlation, Mil further attempted to draw connections be- 
tween the number density evolution of Mg II absorbers n{z) 
and the cosmic star formation history of the universe as 
characterized by the comoving [O II] luminosity density £. 
Given a good agreement between n{z) and i, the authors 
argue that outflows are the mechanism responsible for Mg II 
absorption. 

The Monte Carlo simulation presented in Section 2 also 
allows us to address the agreements in these measurements 
under a simple, generic halo model. It is clear from Figures 
3 & 4 show that both the frequency distribution function of 
Mg II absorbers and the galaxy [O II] luminosity function are 
well reproduced in our mock galaxy and absorber samples 
with no preference in starburst systems. Although under 
the simple halo model Mg II absorbers do not directly probe 
star-forming disks like high column density damped Lya 
absorbers (e.g. Wild ct al. 2007), they trace the halo gas of 
distant galaxies. If more massive galaxies are surrounded by 
more extended Mgll absorbing gas (Equation 7), then n{z) 
is calculated according to 



i (7ri?gas*) 



n{z) 



Ho ^fiM(l-h«)3-FnA 



dL $(L, Z) Kgas fXgas (8) 



where c is the speed of light, $(Lb,2) is the galaxy lumi- 
nosity function, Kgas is the incidence of extended gas that 
is the product of halo gas covering fraction and the fraction 
of galaxies with extended gaseous halos, and (Tgas(i) is the 
cross section of the gaseous halo. As mentioned in § 2.2, the 
mean covering fraction of Mg II absorbing gas is found to be 
high, roughly 100% for absorbers of VKr(2796) ^ 0.3 A at 
0.4 7igas (Chen et al. 2010a). Such high covering fraction 
is consistent with the analysis presented in Kacprzak et al. 
(2008) for a scaling index of /3 = 0.35. 

Combining Equations (4) and (7) yields 



Crgas(I') 



'^gas (tT -Rgas*) 



-^[Oii]* 



-[Oil] 



Kgas (tT Rgas*) 



i[OII]* 

Substituting Equation (9) into Equation (8) leads to 



^loiii- (9) 



n{z) 



e(2)L[oii]* J 



Aoii](2), 



(10) 



e{z)L[ou]* 

where e[z) is defined so that l/e{z) = (c/Ho)(l + 
a) W"m(1 + 2)3 Oa. The same redshift dependent fac- 
tor is defined in Equation (7) of Mil. Equation (10) shows 
that the number density evolution of absorbers naturally 
follows the comoving [Oil] luminosity density, £[ou]{z), if 
extended gaseous halos are a common and generic feature of 
field galaxies. One can therefore apply n{z) for constraining 
the cosmic star formation history, or apply known co-moving 
luminosity density for constraining the fraction of absorbers 
originating in galactic halos (Chen et al. 2000). 

In practice, as illustrated in Mil the SDSS 3" fibers 
define a survey volume of [O II] emitting galaxies. The co- 
moving [Oil] luminosity density, €[0111(2), can be estimated 
according to 

AI/[oiI] 



i;[oii] 



1 



1 



AL 



[Oil] 



AN 



(l + z)'^ Dl{z)dn AN Ar 

Az L[ou] AN 

{l + z)'^ArD\{z)dil Az 



n{z), 



(11) 



[Oil] 



where AI/[oii] is the total [Oil] luminosity observed in a 
coadded QSO spectrum, AA^ is the total number of Mgll 
absorbers over the survey redshift pathlength A z, and L[o ii] 
is the mean [O II] luminosity averaged over the entire galaxy 
(absorber) population. It is clear that Equation (11) resem- 
bles Equation (10) and that £[oii]{z) can be estimated based 
on the product of n{z) and the mean [0 11] luminosity sur- 
face density averaged over all Mgll absorbers. 

We present in Figure 4 the mean [0 11] luminosity sur- 
face density Y^, measured in the 3" SDSS fibers and 

•^^[OII] 

averaged over the number Mgll absorbers in each Wt-(2796) 
bin (solid circles). The mean values can be characterized by 



'[Oil] 



: O3// 



Wr(2796) 



lA 



(12) 



where 03// = 6.6 x 10^^ ergs~^ kpc^ and 63// = 1.1 (solid line 
in Figure 4) . This correlation applies to a 3" survey cylinder 
at 2 = 0.4 — 1.4. Note that because of the extreme bimodal 
distribution of , in which a large fraction of absorbers 

have y^, =0, the median value (Y^, )mod is ex- 
pected to be significantly smaller than the mean X]z,[o nj ' 

While Equation (11) provides a useful tool for estimat- 
ing the co-moving star formation rate density p* based on 
the observations of Mg II absorbers in SDSS QSO spectra (as 
shown cleverly in Mil), we note that the observed strong 
Wr(2796) vs. correlation in SDSS data is shaped 

primarily by a differential fiber loss of the observed [O II] 
flux and therefore unphysical. To illustrate the fiber selec- 
tion bias in driving the apparent correlation in Figure 4, we 
repeat the Monte Carlo analysis described in § 2 for 15" di- 
ameter fibers. This is five times the aperture size of SDSS 
spectra. At 2 « 1, such a large aperture covers an area of 
roughly 40 kpc projected radius from QSO lines of sight, 
sufficient to cover the majority of the light from galaxies pro- 
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Figure 5. The same as Figure 4, but assuming 15" fibers in the 
Monte Carlo simulation described in § 2. Both )mcd (fed 

stars) and (solid blue circles) show a lack of correlation 

with iyr(2796). The bottom panel shows that with 5x larger 
fibers, differential aperture loss has minimal impact on the survey 
of Mg II selected galaxies along QSO lines of sight. 



ducing strong Mgll absorbers (Wr(2796) ^ 1 A). Figure 5 
shows the distribution of Wrf2796) and Y^, measured 

over 15" diameter fibers for the mock Mgll absorber sam- 
ple. The mean [Oil] fiux X/L[o ] ^^^^^ ^ similar slope as 
the median [Oil] flux )mGd within the 15" diameter 

aperture and the slope is consistent with 615// = (sohd line 
in Figure 5). 

In summary, we have demonstrated that (1) because 
of the observed Wr{2796) vs. p anti-correlation (e.g. Chen 
et al. 2010a), galaxies at larger projected distances pro- 
duce on average weaker Mgll absorbers and weaker (or 
zero) i[oii] in SDSS QSO spectra and that (2) because 

of the extreme bimodal distribution of Y^, in which 

^-£-[0 11] 

a large fraction of absorbers have Y^, =0, the me- 

'^-'-[OII] 

dian value (X^^^^ j)med can be significantly smaller than 



the mean ^ 



[O II] 



parent Wr{2796)-{J2L^^^^^/^cc 



Together these effects strenthen the ap- 
j)mcd correlation of Mil. 
Consequently, empirical correlations between star- 
forming properties of galaxies and statistical properties of 
QSO absorbers do not provide the unambiguous evidence 
necessary to discriminate between infalling clouds and out- 
flows as the mechanism for producing Mg II absorbing clouds 
at p > 50 kpc of a galaxy. While outflows are a natural 
product of starbursts, gas accretion provides the fuels for 
star formation in galaxies. As discussed by previous authors, 
mass is likely a more fundamental factor that determines the 
properties of gaseous halos around galaxies (e.g. Steidel et 
al. 1994; Ledoux et al. 2006; Chen et al. 2010b) and more 
massive galaxies can support more extended gaseous halos. 
Given that models that do not require outflows can also 
reproduce the empirical correlations between absorber and 



star formation properties, we caution drawing conclusions 
in favor of outflows based on simple correlations between 
absorber and star formation properties. 
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